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Abstract 



We present a generic solution to the fundamental problem of how to con- 
nect two points in a plane by a smooth curve that goes through these points 
with a given slope. The smoothness of any curve depends both on its cur- 
vature and its length. The smoothest curves correspond to a particular 
compromise between minimal curvature and minimal length. They can be 
described by a class of functions that satisfy certain boundary conditions 
and minimize a weight functional. The value of this functional is given es- 
sentially by the average of the curvature raised to some power v times the 
length of the curve. The parameter v determines the importance of minimal 
curvature with respect to minimal length. In order to find the functions that 
obtain the minimal weight, we use extensively notions that are well-known 
in classical mechanics. The minimization of the weight functional via the 
Euler-Lagrange formalism leads to a highly non-trivial differential equation. 
Using the symmetries of the problem it is possible to find conserved quanti- 
ties, that help to simplify the problem to a level where the solution functions 
can be written in a closed form for any given v. Applying the appropriate 
coordinate transformation to these solutions allows to adjust them to all 
possible boundary conditions. 
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I. INTRODUCTION 



A. The problem 

Consider two points in a plane, each associated with a ray pointing in some direction. 
The problem we address in this paper is how to connect these points in a "smooth" way. 
That is we are looking for a curve that satisfies the following three conditions: 

1. The curve goes through both points. 

2. The associated ray at each point is a tangent to the curve. 

3. Between the two points the curve follows some optimal path, which is a compromise 
between minimal curvature and minimal length. 

The first two requirements constitute the boundary conditions of the problem. The third 
condition defines what we mean by "smooth" . 

B. Motivation and Outline 

A solution to the problem outlined above has many obvious applications. For example it 
could be used to determine the ideal shape of a road to be built between two points where its 
direction is predetermined (e.g. by two bridges). In fact any interpolation problem that has 
been reduced to the task of connecting a sample of data points where the slope is fixed can 
be solved using the elementary solution of connecting smoothly two points. If the original 
sample only consists of points, the respective slopes can be determined, for example, by 
taking the slope of the line through the neighboring points, or by some different, more 
sophisticated method. Evidently the quality of interpolation curves is crucial to all fields 
that deal with numerical data from applied sciences to economics. It is interesting to note 
that while there exist many interpolation schemes, most of them rely on simple functions, 
like polynomials, that in general do not give the best interpolation. The aim of this work 
is to establish first a criterion for the quality of the interpolation and then to investigate 
quantitatively which are the optimal interpolating curves. 

As we shall see in the following, the problem outlined above translates into a well-defined 
mathematical exercise, which is interesting by itself. In fact we feel that our solution to the 
problem, relying heavily on notions well known in physics, has also a pedagogical value. 
It serves as an example of some fundamental physics principles in the context of a very 
intuitive and visual, yet non-trivial problem. 

Since the problem is defined as a geometrical task, its mathematical formulation ob- 
viously does not depend on time. Instead we take one space-direction (x) as the variable 
of integration, and the other space-direction (y) to describe the solution function. While 
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this description is not manifestly invariant under rotations and translations in the x — y 
plane, it allows us to define an action functional that assigns a weight to any function (that 
complies with the boundary conditions) according to its "smoothness". It is given essen- 
tially by the average of the inverse curvature radius, to some power u, times the length of 
the curve. This action is minimized by a class of functions, each describing the smoothest 
curve for a particular choice of the parameter v that specifies the relative importance of 
minimal curvature with respect to minimal length. These functions are the solutions of the 
Euler- Lagrange equation, which turns out to be a complicated, non-linear third order differ- 
ential equation. Solving this equation can be facilitated immensely by taking advantage of 
the symmetries of the problem. We compute the linear and angular momenta that follow 
from Noether's theorem due to the translational and rotational invariance. The equation 
of motion is also scale invariant, but there exists no conserved charge corresponding to 
scaling, since the Lagrangian and other dimensionful variables do change under scaling 
transformations. The problem is solved explicitly for a particular choice of the momenta. 
Applying the appropriate coordinate transformation to this solution allows to adjust it to 
all possible boundary conditions. 



II. FORMALISM 



In order to translate the problem outlined in Section |I A| into a mathematical one we 
have to introduce some notations: Let us give "names" to the two points: We shall refer 
to them as the initial point Pi and the final point Pf. Which one is which is arbitrary, but 
when introducing coordinates 

Pi : (xi,yi) and P f : (x f ,y f ) , (1) 

we demand that X{ < Xf. The coordinates (x,y) are a set of two real numbers specifying 
any point in the plane. We use simple Euclidean geometry. The associated rays at Pi and 
Pf are described by their inclination angles, cti and a/, with respect to the vector pointing 
from Pi to P f (c.f. Fig. |1|). 

We call the function y(x) that satisfies the three conditions in Section [IA| the "smooth 
connection function" (SCF). Our task is to determine this function. The boundary condi- 
tions defined in Section [TS| are: 

y(xi) = yi and y(x f ) = y f , (2) 

y'(xi) = tan(«j + ao) and y'{%f) — tan(a/ + «o) , (3) 

where the prime denotes a derivative with respect to x and a is the angle between the 
positive x-axis and the vector pointing from Pi to Pf. 

In order to include the third and crucial condition in Section ITS] we need to incorporate 



both the curvature and the length of the curve. The local curvature at a given point on 
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the curve is defined as the inverse radius (1/r) of the circle that coincides with the curve in 
the infinitesimal vicinity of the point. It is given in terms of the first and second derivative 
of the function y(x) by: 

1 _ (4 ) 



r(x) [l + y'(x) 2 } 3 / 2 ' 

To show this relation it is enough to verify that the curves that describes a circle, i.e. 



y(x) = ±yV 2 — (x — xo) 2 + yo, where r is the radius and (xo, yo) denotes the center of the 
circle, give r(x) = r. The length of an infinitesimal piece of the curve at x is given by 



ds(x) = Jl + y'(x) 2 dx. (5) 
We define now the following functional: 

S[y(x)} = [*' ds{x) = H |y/,(g)r 3y _ 1 dx. (6) 



r(x) 



We shall refer to S[y(x)} as the action. For a given function y(x) it returns the weighted 
sum of ds(x) where x goes from Xi to Xf. The weight at position x is given by the local 
curvature l/r(x), defined in eq. (|J), raised to some power v. A priori this parameter is 
an arbitrary real number that determines how to choose the compromise between minimal 
curvature and minimal length. For example, for v = the curvature does not play a role 
at all and S[y(x)] u= o is just the length of the curve. For \v\ — > oo the situation is exactly 
the opposite, since S[y(x)]\ ll \^ 00 — J (l/r) u dx only depends on the curvature in this case. A 
special situation arises for v — 1, where S[y(x)] = J^/(ds/r) = f* } da = af — ai. Since the 
action gives just the difference between the inclination angles at the boundary, independent 
of the particular choice of y(x), it is impossible to determine the optimal path for v = 1 
and we exclude this case from our subsequent discussion. 

In order to find the optimal curve we have to minimize the action under the boundary 
conditions in eq. (0) and eq. (|3|). Before we continue with our analysis we would like to give 
a concrete example illustrating the relevance of the action functional. Consider a spacecraft 
that has been launched to explore some distant planets and which is flying with constant 
velocity v through the interstellar space. When approaching a certain planet, one has to 
adjust carefully the trajectory of the spaceship, say to enter the atmosphere of the planet 
at a precise angle or to use its gravitational field to accelerate the spacecraft to its next 
destination. Suppose that for this purpose the spaceship has thrusters that exert a force 
F perpendicular to the direction of motion. Then at any given time t between the ignition 
(at tj) and the end of the manoeuvre (at tf) the instantaneous curvature radius r(t) is 
inversely proportional to F. Let us assume that the fuel consumption per unit time, W, is 
governed by some potential law, i.e. W oc F u oc (1/r)", where v is an empirical parameter. 
Then the total fuel consumption, C = J t f Wdt, is proportional to the action in eq. (H), 
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since dt = ds/v. Obviously minimizing C is crucial. Once we have found the solution for 
the trajectory we can obtain the curvature radius as a function of t which indicates how 
much fuel should be burned at each point. 

The integrand in eq. (^J), usually referred to as the Lagrangian, 

c.{y',v")= ly " {x)l \^ (7) 

only depends on y'(x) and y"(x), while there is no explicit dependence on y(x) and x. Thus 
y and x enter the action only via the boundary conditions. From eq. (0) it follows 

that 

f x f 

Vf — Vi = V( x f) - V( x i) = / V ( x ) = constant . (8) 

JXi 

This constraint can be absorbed into the action functional by introducing a Lagrange 
multiplier A, resulting in a new action 

S[y'{x)\ = S[y'{x)\ + A P y\x) = P [I + A y\x)\ dx . (9) 

Since there is no reference to y(x) anymore, but only to its first and second derivative, we 
can change variables, 

q(x) = y'(x) and q'{%) — y"{ x ) , (10) 
and write the Lagrangian corresponding to the action <S[g(x)] in eq. (0) as 

C(q, q') EE C(q, q') + A q(x) = + A q . (11) 



1 + 9 



2^! 
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Now, to minimize S[y(x)] under the boundary condition in eq. (^|) is equivalent to the 
minimization of <S[g(x)]. A necessary condition for <S[g(x)] to be extremal is that its first 
functional derivative SS[q(x)]/Sq(x) vanishes. Since <S[g(x)] only depends on q and q', this 
is a well-known problem. Its solution is given by the Euler- Lagrange equation: 

8C d9C 0. (12) 



Computing 



dq dx dq' 
dC (l-3u)\q'\ u q 



dq (! + q 2 
dC va\q'\ v - 1 



3g+l 



Zv-l 1 



dq' (! + q 2 
d dC v{y- l)\q'\ v ~ 2 q" v{\ - 3v)\q'\ u q 



+ A, (13) 

(14) 
(15) 



dx dq' ( i + g 2)— ( 1 + g 2^ — 



where a = sign(g'), we obtain 

x (f-l)l? 



+ g 2 )?" + (l-3,)( 9 ') 2 d . (16) 



■t\ 2 



(1 + 9 

The above equation still depends on the parameter A. We can eliminate this dependence by 
differentiating eq. (|T6|) with respect to x yielding the following equation of motion (EOM) 



= {V 1)|g 3 l + 3~ {(1 - 31/) [(1 - 3^ 2 )(g') 4 + 2vq{\ + q 2 ){q'fq 
{1 + q 2 )— L L 



/\u-3 

'\2ll 



+ 

v{l + q>Y[{p-2){q"Y + q'q">]}. (17) 



For all functions q(x) that solve the above EOM the action 5[g(x)] in eq. (|9|) is stationary, 
i.e. 5S[q(x)]/5q(x) = 0. However, in order to minimize S[q(x)} also the so-called Legendre 
condition, 



d 2 C u(u~l)\q'^- 2 
0< !U 



WY (1 + q 2 ] 



3v-l J 
2 



has to be satisfied. This condition arises from the second functional variation of the action 
with respect to q(x). Note that unlike for the minimization of usual functions it is only a 
necessary condition for a minimum of 5[g(x)]. This is essentially because a small variation 
in q does not necessarily imply a small variation in q' . However in our case the sign of 
d 2 C/(dq') 2 only depends on the prefactor z/(z/ — 1), but is independent of q and q' even 
if they are taken as independent variables. Such a situation is referred to as a "regular 
problem" and it implies that the solutions of the EOM indeed give rise to a local minimum 
of the action if v < or v > 1. Conversely for < v < 1 one obtains a local maximum. 
Thus we rule out any < v < 1 from our analysis and we shall only consider v < and 
v > 1 from now on. 

Even though we have managed to translate the problem at hand into a differential 
equation, solving this equation analytically seems like a formidable task given that eq. ( |17D 
is third order in q and non-linear. However, the integrated form in eq. ( |T6"D gives us a hint 
of how to proceed. The fact that A is a constant is related to the symmetries inherent 
to the problem. Studying these symmetries systematically one can simplify the problem 
significantly and find its solutions. 

III. SYMMETRIES 

Let us examine more carefully the problem with respect to its intrinsic symmetries. 
The important point to realize is that when introducing the coordinates in eq. (|l|) we have 
made an explicit choice of 
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• where to define the origin, 



• how to orientate the axes and 

• what units of length to use. 

These choices are arbitrary and once we have found a solution to our problem we can 
redefine our coordinate system. Assume that y(x) solves the EOM in eq. ([17|) in some 
coordinate system E and that y(x) is obtained from y(x) by a transformation to a different 
coordinate system E , 

(A fS\ R fx\fz \ ^ (19) 



v \y \y \vq 



where 



y sin cos J 

describes a rotation by an angle 8, and Xq and yo parameterize a translation in the Mis- 
direction and in the y-direction, respectively. The action in eq. (^) is defined in terms 
the curvature radius r [c.f. eq. (f|)] and the length element ds in eq. (|5p. Both quantities 
are manifestly invariant under rotations and translations. Therefore it is clear that the 
coordinate transformations in eq. (|T9| ) do not change the action nor the EOM derived from 
it. However the boundary conditions obviously have to change under E — ► E. 
Let us consider the infinitesimal change of the coordinate y at a given x, 

$y{x) = y(x) - y(x) , (21) 

for the three types of coordinate transformations. For a translation in the y-direction 

yo = £y<l, ^o = 0, 9 = 0, (22) 

the answer is trivial. The change in y is simply 

5y(x,e y ) = e yj (23) 

implying that the changes in all the derivatives of y(x) vanish: 

Sq(x,e y ) = 0, 

5q'(x,e y ) = 0, (24) 
5q"(x, €y) = 0. 

However, a translation in the x-direction 
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2/0 = 0, x = e x < 1 , 
is a bit more tricky. First, note that 

y(x) = y(x) . 



6 = 



(25) 



(26) 



This expresses merely the fact that the coordinate transformation in eq. (pi]) leaves the 
functional behavior of y(x) invariant. It does not make a difference whether we refer to 
the function by y(x) in the coordinate system £ or y(x) in the system E. Then a Taylor 
expansion of y(x) around x gives 



y{x) 



y{x) + Tx 



(27) 



implying that 



5y(x,e a 



dy_ 

dx 



(28) 



From eq. fl28|) it follows that the corresponding changes in the derivatives of y(x) are given 
by 



5q(x, e x ) = -q'(x) e x , 
5q'(x, e x ) = -q"{x) e x 
8q"{x,e x ) = -q"'{x)e x 



(29) 



Finally, an infinitesimal rotation 

2/o = 



x = , 9 = e e 



(30) 



involves both a change in x by — yeg and in y by xeg. Together this results into a change 
of y{x) by: 



5y(x,€ ) = [x + 



dy_ 

dx 



(31) 



From eq. ( |3~lD it follows that the changes in q, q' and q" for an infinitesimal rotation are 
given by 



5q(x,e e ) = (1 + q 2 + yq') e e , 
5q'(x,e e ) = (3qq +yq") e e , 
5q"{x^e) = {Z{q l ) 2 + ±qq" + yq l ")e e . 



(32) 



Let us now investigate the effect of the coordinate transformations in eq. (|19|) on the 
Lagrangian C in eq. (|TT|). In general an infinitesimal transformation parameterized by 
e< 1 will induce a change 
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8C = £(q, q') - C(q, q') 



(33) 



Provided that C has no explicit dependence on x, a Taylor expansion of C(q, <f) around q 
gives 



8C = —8q + —8q' (34) 




?£ - — f?£\ s q (35) 
dq dx dq' J 

The second term in eq. (^) vanishes for q(x) that satisfy the equation of motion (|T^) . 

Any transformation that changes the action iS[g(a;)] in eq. (^) at most by constant is 
called a symmetry transformation. Such a transformation does not change the extremal 
condition for the action and therefore leaves the EOM invariant. (Note however that not all 
transformations that leave the EOM unchanged are symmetry transformations according 
to our definition.) Then, the corresponding Lagrangian could change only by a total 
derivative: 

dF 

5C = e— , (36) 
dx 

where F(q, q') is a function of q and q' . Equating the two results for 8C we find that 

dF_ = d_ fdC \ 
dx dx \dq' J 

We stress that the expression on the left-hand side is correct for a symmetry transformation 
of any function y(x), while the right-hand side is only valid for a transformation of y{x) 
that solves the EOM. From eq. ([37]) it follows that the conserved charge defined as 

is a constant with respect to x, i.e. dQ/dx = 0. Consequently it has one specific value for 
all points on a particular solution of the EOM. Obviously Q can only be defined up to a 
constant. The above argument, that every continuous symmetry transformation implies a 
conserved charge is known as the Noether theorem. 



To be explicit let us apply it to a translation in the x direction. Using eq. (p9|) it follows 
that 



8C(x,e x ) = -^-q+-q je x 



(39) 



-7-e x . (40) 
dx 
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This means that (up to a constant) we can identify F x = —C. Then it follows that the 
conserved charge corresponding to translations in the x direction is given by: 



dC 



dq 



7? 



C 

i)k 



— - Xq. 



(41) 



We call P x the conserved momentum in the x direction. In fact it is nothing more than 
the Legendre transformation of the Lagrangian with respect to q' . (If the variable of 
integration in eq. (§) had been the time t rather the coordinate x then the equivalent 
Legendre transformation of the Lagrangian with respect to dq/dt is called the Hamiltonian 
and the charge related to time invariance is the energy.) 

One might guess that a similar argument for translations in the y direction should 
give another conserved charge P y , which is the conserved momentum in the y direction. 
However a change in the variable y(x) as in eq. ( ^231) has no effect on q(x) and q'(x), see 
eq. fl24"|), implying trivially that 



5C(x,e y ) 



0. 



(42) 



Thus Noether's theorem does not help in this case to derive the conserved charge. However 
we have already encountered another conserved quantity, which could serve as a candidate 
for P y . In order to absorb the boundary condition into the action we introduced the 
Lagrange multiplier A, which intuitively is related to changes in y, c.f. eq. (H). Eq. ( |i~6l) 
states that A equals to a complicated function of q, q' and q" for all x. Therefore this 
function is a constant of motion. To prove that indeed P y = A is nontrivial and we will 
show this after discussing the conserved charge related to the rotations in eq. flT5|). 

Using eqs. flilf ) and ( |i~4"D and the changes in q and q' under rotations according to 
eq. fl32|) the change in the Lagrangian under an infinitesimal rotation is 



5C(x,e e ) 



[i-3v)\g7q , y 

3U+1 r A 

(1 + g 2 ) — 



[l + q 2 + yq') 



v\q 



l\v-\ 



[1 + q 



3v-l 



(3qq' + yq")\e e . (43) 



It is not difficult to check that this can be rewritten as a total derivative, 



SC(x,e g ) 



d_ 

dx 



y\Q 



— - + X(x + yq) 



(44) 



;i + g2) . 

Thus for infinitesimal rotations the function F can be identified to be (up to a constant) 

F e = yC + xX . (45) 

Then the charge corresponding to the rotation symmetry, which is called the total angular 
momentum, is given by 
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J 



dC 5q(x,e e ) 



dq' 
yC + xX — 



v a\q 



l\v-\ 



(I 



3v- 
2\ — 



- (1 + q 2 + yq') 



(xP y -yP x )+S. 



(46) 

(47) 
(48) 



In the last step we have separated the total angular momentum into the orbital contribution 



L — xPy yPx 



and the remaining term 



S = - 



v a\q 



l\V-\ 



v a 



(1 + q^ Hx)]»-i> 



(49) 



(50) 



which we call the spin or intrinsic angular momentum. The orbital angular momentum 
would be the sole contribution if y(x) is a straight line (as is the case for v = or if q' = 0). 
A non-vanishing spin S arises for all the other curves due to their curvature. Note that it 
is sensitive to the sign a of q' . 

We would like to come back now to our claim that P y coincides with A defined in 



eq. fllB]) . To this end let us compute the change in Pj (i 
general 



x, y) induced by rotations. In 



8PAx, e e ) 



BP- 

f< 1 + ^ 



(51) 



where we used the results in eq. fl32[) and the fact that the terms proportional to y add up 
to the total derivative of P,-: 



dq dq' dq" 



dPi 
ax 



(52) 



Assuming that indeed P y = A, a somewhat tedious but straight-forward calculation of 5P y 
according to eq. fl5Tp yields 

{y-l)\u 



l\u-2 



5Py(x,e e ) 



[l + q 



3i/+l 
2^1 — 



1 + 3uq')(q') 2 - uq(l + q l )q" e e = P x e e . 



2\ _// 



(53) 



Similarly using eq. (|5lD one shows that SP x (x, eg) = —P y ee- From the infinitesimal transfor- 
mations it is clear that P = (P x , P y ) transforms like a vector under rotations, i.e. P = R P, 
where R is the rotation matrix defined in eq. fl2"0"|). It follows that the constant A in eq. ( |I~6"D 
can indeed be identified with the conserved momentum in the ^/-direction. 

Using the above results it is easy to compute the change in the orbital angular momen- 
tum L under infinitesimal rotations: 
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6L(x, eg) = x ■ SP y (x, e e ) - [Sy ■ P x (x, eg) + y ■ SP x (x, e e )] 

= -y(qP x - P y )e e (54) 

Note that this is the change at a fixed x, so of course there is no variation with respect 
to x. The change in the spin S under rotations is given by 

dS OS 
5S(x,e e ) = —5q(x,e e ) + —5q'(x,eg) 



v{v - 1 



~3q(q>) 2 - (1 + q 2 )q" 



3i>-l 

,2\—r~ 



eg 



(l + q 2 

= y(qP x - P y )e e . (55) 

It follows that the changes in L and S exactly cancel each other such that the total angular 
momentum J does not change under rotations, i.e.: 

8J(x, eg) = 5L(x, eg) + 5S(x, eg) = . (56) 

For completeness let us also compute the changes of the three conserved charges under 
infinitesimal translations. Since the momenta P x and P y do not depend explicitly on y(x) 
they are trivially invariant under translations in the ^/-direction by Sy(x, e x ), due to eq. (p4j). 
Infinitesimal changes in the x-directions induce Sy(x, e x ) as given in eq. ( p8] ) implying that 
the change of the n th derivative of y(x) is given by ?/ n )(x, e x ) = —e x y <yn+l \x), c.f. eq. (^). 
Then the changes in the momenta are proportional to their total derivative, which vanishes, 
i.e. 

(dPi , dPi „ dPi ,„\ dPi 
5PAx, e x ) = - —^q' + — -q" + — -q'" e x = -^e x = . 57 

Finally the total angular momentum J in eq. ( |4lj|) does depend explicitly on y(x) and 
consequently changes under translations in the ^-direction by 

8J(x, e y ) = -P x e y . (58) 

For the change in J due to an infinitesimal translation in the x-direction we obtain an 
expression similar to the one in eq. (|57|): 



. T/ . (dJ dJ . dJ .. dJ ,.A / dJ d.J\ _ . 

* e ■> = " ( ai q+ V + ®t q + W q ) t '--\T x -Tr x ) t ' = ^ ■ (59) 

where we used the fact that dJ/dx = and that J depends explicitly on x. We note that 
the results in eqs. (59) and (|58|) are also correct for finite translations, since the changes 
do not depend on J. 
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IV. SCALING AND DIMENSIONAL ANALYSIS 



Consider a scaling transformation that changes x and y by a fraction a of their original 
value, i.e. 



x 

y 



X 

V 



X 

y. 



+ a 



x 

y. 



(60) 



The infinitesimal change in y(x) at a fixed x under such a transformation for a = e„ <C 1 
is given by 



dy_ 

dx 



5y(x,e a ) = \y 
Consequently the derivatives of y(x) change by 



x 



(61) 



Sq{x,e a ) 
Sq'(x,e a ) 



8q 



iq'x) e a , 

-{q"x + q) e a , 
-{q {n+l) x + nq {n) )t a . 



(62) 



Therefore it follows that some arbitrary function A of y(x) and its derivatives changes 
under an infinitesimal scaling transformation by 



SA(x,e c 



' (dA 


dA\ 


\ dx 


dx J 



OA 



OA 



X dy y + ^dq^ 



nq 



(«) 



(63) 



Using this formula it is easy to show that the conserved charges P x , P y and J transform 
as follows under the scaling transformation: 



SPi(x,e c 
5J(x,e c 



x,y) 



(64) 
(65) 



Note that for each charge the infinitesimal change is proportional to its original value. 
This implies that the finite scaling transformations are given by Pi Pi - exp(— va) and 
J — > J ■ exp[(l — u)a}. 

The interpretation of the proportionality factor follows from the following argument. 
Any quantity A can be written as a product of a dimensionless number and the dimension 
[A]. Since the problem we discuss is purely geometrical we only have a fundamental 
length scale Iq. So [A] = Zg can be written as some power of this scale. Now the scaling 
transformation in eq. (^) can be viewed as a change of the length scale Iq by some fraction 
51q = e a /o- Then to first order the corresponding change in A is given by 



5 A = —5l 
dl 



[i A e a . 



(66) 
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Thus the proportionality factor is nothing more than the dimension of the quantity and 
we find that [P x ] = [P y ] = Iq v and [J] = l ~ v - It is reassuring that these results follow 
also from dimensional analysis of the definitions of the conserved charges by noting that 
[x] = [y] = l and [q^] = Iq u . The dimension of the action [S] = ll~ u allows us to 
understand a posteriori why the regime < v < 1 had to be excluded from our analysis. If 
v = 1 the action is invariant under scaling transformations. In particular, we can shrink or 
magnify sections of any possible solution and thereby transform it to any arbitrary shape 
without changing its action. This explains why for v = 1 the action only depends on 
Act = a.f — ati, as mentioned after eq. (|6]). The physics jargon is to say that the action 
becomes "soft" when v approaches unity and it is "critical" at v — 1. Now if < v < 1 one 
can always find a "trivial solution" which is defined as follows: just follow the rays at P, 
and Pf to the point where they intersect and bend the curve in the infinitesimal vicinity of 
the intersection point by Act. For this curve the action vanishes, since any length element 
ds(x) of the straight part of the curve, where r(x) = oo, does not contribute to the action 
as long as v > 0. Moreover the (infinitesimal) part of the curve that is bended also does not 
affect the action, because from the scaling property of the action we know that it decreases 
when shrinking the unit length Iq as long as v < 1 . Consequently continuously scaling down 
the region where the bending takes place, we achieve a zero, and hence minimal action. 
For v < the fact that any finite straight piece of the curve gives an infinite contribution 
to the action precludes this solution and for v > 1 it is not viable since, due to the inverse 
scaling, the action blows up if the curve is bended strongly on a section of small length. 

The fact that the momenta and the action have somewhat unusual (and v dependent) 
dimensions could easily be remedied by multiplying the Lagrangian in eq. ([11]) by Iq~ . 
Then the action and the angular momenta would be dimensionless and the linear momenta 
would have dimensions of Iq 1 . 

Using eq. ( |63| ) or just applying dimensional arguments it follows that a scaling trans- 
formation on the right-hand side of the EOM in eq. ( |T7D results in a multiplication by a 
factor [1 + (1 — u)e a ]. However since the left-hand side is zero it is clear that the EOM is 
unchanged in the new coordinate system. The important point to note is that even though 
the EOM is invariant under scaling, the transformation in eq. (|60|) is not a symmetry trans- 
formation. The reason is that the change in the Lagrangian under an infinitesimal scaling 
transformation, 



cannot be written as a total derivative for v ^ 1. As a consequence there is no conserved 
charge related to scaling. 




(67) 
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V. THE SOLUTION 



The coordinate transformations discussed in the previous section are very useful for 
actually solving our problem: If we manage to find the SCF in some convenient coordinate 
system £, we can use translations, rotations and scaling to fit the particular solution to any 
boundary conditions given in some other coordinate system E. Let us consider a particular 
solution for which P x = J = 0. Then from eqs. fl4T|) and (|48|) it follows that 

xPy - ^^al • (69) 

2 



From eq. (B9|) we get 



M = (f?£)~ ( l + (70 ) 

provided that axP y / v is non-negative. This requirement implies that when x changes sign 
also q' has to change its sign a. Plugging the result for \q'\ into eq. (|68|) and solving for q 2 
we find 

✓-ii i — 

2 C\X\—^ 

1 =~ ^, ,_2^ ' ( 71 ) 

1 — C \x\ "-i 

where C = \y — 1) |i / | T=J '|Pj / |' ;= i can be set to unity by choosing P y = ±|vp/|f — lp - . 
Because g can be positive or negative there are two solutions 



y(a?)± = ± 



+ y±(x ) , (72) 



x 



2^ „-! 



where y±(xo) are the respective integration constants. Since either v < or v > 1 the 
solution in eq. ( [72] ) is defined only in the interval —1 < |x| < 1. 
In the limit where \u\ -> oo we can solve eq. (|72|) analytically: 



lim y(x) ± = ±[ X J ^— 2 dx = ± (Vl-zg - Vl - x 2 ) + y ± (x ) . (73) 

We see that in this case the SCF describes a segment of a circle, which has a constant 
curvature radius and therefore presents the best solution if we only care about minimal 
curvature along the curve. For finite values of v also the length of the curve plays a role. 
The integral can be expressed in terms of hypergeometric functions F(a, b; c; x), i.e. 
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V(x)± = Txlm [F{^, i; ^; + y±(z ) . (74) 

We show the functions y(x) + (solid) and (dashed) for various z/ in Fig. [| The two 

branches are monotonic and we have chosen the constants of integration y±(xo) such that 
they can be obtained from each other by a reflection with respect to the y-axis. (If one 
chooses the constants of integration such that all curves go through (1,0) the two branches 
are related by a reflection with respect to the x-axis.) Each SCF changes the sign of its 
curvature at x — 0. For large \v\ the SCF is very close to the arc of a circle. The curves 
for v > 1 are below the curve for \v\ — > oo and they become "flatter" and thus shorter 
for smaller values of v. The curves for positive v approach x = with a vanishing slope 
and their curvature radius diverges at x — 0. The closer v is to unity the sooner the curve 
approaches the value y(0) when x — > (which can be understood from the scaling behavior 
of the action c.f. section [TV]) . The curves for v < all reside above the curve corresponding 
to \v\ — > oo. For these curves both the slope q(x) and the curvature radius r(x) vanish at 
x = 0. 

The standard smooth connection functions shown in Fig. |2] are the fundamental so- 
lutions to the boundary problem we want to solve. Any specific solution consists of a 
segment of a standard smooth connection function that can be viewed as a template which 
may be rotated, translated and scaled in order to fit the boundary conditions. Before we 
continue to describe in detail how this can be done, it is useful to extend the SCF beyond 
the interval they are defined on. 

To this end we note that from eq. ( [70] ) it follows that the curvature radius, defined in 
eq. (fD, is given by 

r{x)± = \x i-" • . (75) 

v 

It has the same value for the two branches of the SCF in eq. (|72|) . Therefore, even though 
q(x) diverges at \x\ = 1 the curvature radius has a well-defined limit for \x\ — > 1, 

lim r(x)± = . (76) 

Thus it is natural to connect the two solutions y(x)± to one single curve and for the 
following we shall fix the constant of integration such that \im\ x \- t iy(x) = for all curves. 
The problem is that these curves are not single- valued. In order to express them as a single 
function we simply exchange the coordinates x and y. The resulting curves are shown in 
Fig. |||. All curves have been rescaled such that they are defined in the interval —1 < x < 1. 
The interval —1/2 < x < 1/2 corresponds to the SCFs of Fig. 0. The continuations beyond 
x = ±1/2 use a segment of the other branch, which is attached to the centerpiece such 
that the curvature radius is continuous at x = ±1/2. 

Now that we have these "extended standard smooth connection functions" of a given 
v > 1 for a particular set of values for the conserved charges, namely 
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J = 



and 



1 


V 


V 




v - 


- 1 


u-l 



(77) 



it is not difficult to obtain a specific SCF for any given boundary conditions in eqs. @ 
and (fj). The basic idea is to find first two points Pi and Pf on a "standard SCF", where the 
slopes correspond to the required angles ctj and atf, and then to apply a set of coordinate 
transformations to the curve in order to match the boundary conditions. The first step 
implies that we have to check whether it is possible to find positions Xi and Xf somewhere 
on the standard SCF such that 



q(xi) = tan(«j + <5o) and <?( x /) = tan(a/ + So) 



where 



tan diQ = t (xi,Xf) 



Xf 



x ; 



(7? 



(79) 



Using that tan(a + (3) — (tan a + tan f3) / (1 
as two coupled equations 

tan(ttj) = ti(xi,Xf) : 

tan(«/) = tf(xi,Xf) 



tan a tan (3) we can rewrite these conditions 



q(xj) ~ t (xi,x f ) 
1 + q(xi) t (xi,x f ) ' 
qjxf) - t (xi,x f ) 
l + q(xf)t (xi,x f ) ' 



JO) 
51) 



This set of equations can be solved numerically. For example one can apply Newton's 
method and use the iteration scheme 



~(n+l) 



X 



f 



'x {n) 



dU 


% 


dxi 


dxi 


dU 




dx f 


dxi 



-1 



(Xij) — tan(a,: 



- ( n) \ - tan(a/) 



12) 



Of course such an iterative procedure will only converge provided that for a given v one 
can indeed find two points on the standard SCF that satisfy eq. fl78p. The important 
observation is that using the extended standard SCFs (as shown in Fig. it is possible to 
find a suitable segment of the curves for any given pair of (aj, atf). Thus we will use these 
curves in the following. 

In general the endpoints Pj and Pf of the fitting segments do not satisfy the boundary 
conditions in eqs. (§) and (|3|). However, it is important to realize that neither translations, 
nor rotations, nor scaling transformations change the angles o.i and a/. This is because 
they are defined relative to the line through Pj and Pf. The only variable in eq. ( |7S| ) that 
does change under coordinate transformations is the angle <5o that defines the direction of 
this line with respect to the x-axis. Thus we can apply a rotation in order to match a with 
any given value in eq. (H) for olq. Since a® is defined as the ratio between the differences 
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of the coordinates of Pi and Pj [see eq. (ff9[)1 it is invariant under translations and scaling 
transformations. This enables us to satisfy also the boundary conditions in eq. (0) by first 
scaling the (rotated) SCF y(x) such that the distance between Pi and Pf coincides with the 
distance between Pj and Pf and then translating the resulting curve such that it connects 
these points. The SCF obtained like this satisfies both eq. (fj) and eq. (H). In Fig. | we 
show a specific set smooth connection functions, each corresponding to a different v, that 
all satisfy the same boundary conditions. Various other examples of smooth connection 
functions, that all have been obtained by the numerical recipe described above, are shown 
in Fig. [5] Each plot corresponds to a particular choice of («j, a/) and shows the behavior 
of the solution for five different values of v. 

Finally we would like to compare the SCF discussed in this paper with other "conven- 
tional" interpolation functions. In particular for practical purposes it is important to know 
how much better the optimal path (i.e. the SCF) is with respect to some approximation. 
To this end it is useful to have a closed expression for the action in eq. (pj), which deter- 



mines the "smoothness" of any curve y(x). Using eq. (II]) its integrand can be written in 
terms of the linear momenta as 

If y(x) solves the EOM, then P x and P y are constants and C can easily be integrated 
resulting in 

§ = P ^ x f - x i) + P v(Vf ~ Vi) ( g4 ) 

Note that S is manifestly invariant under rotations and translations. If y(x) does not 
solve the EOM, the integration in eq. (Q) in general has to be done numerically. Standard 
interpolation curves are often given in parametric form, i.e. as {x(t),y(t)}, where t varies 
within a given interval [ti,tf] along the curve. In this case the action is given by 

~ ftf \xy-xy\ v 

S p = / 777 T^EX dt > ( 85 ) 

J u {x z + y z ) 2 

where the dot denotes a derivative with respect to t. We note that in principle one could 
also minimize this action in order to determine the smooth connection functions. However 
in this case the Euler-Lagrange formalism results in two coupled non-trivial differential 
equations, which appear to be even harder to solve than the EOM we got in eq. (|17|). 
Using eq. (^) we have computed the weight functional for cubic Bezier curves. In general 
for small v the SCF can be well-approximated by an appropriate Bezier curve and the 
respective values for the action only differ at most by a few percent. However, for larger v 
it is increasingly difficult to match the SCF by a Bezier curve and the difference between 
the actions becomes significant. This is because for v ^> 1 the SCF approaches an arc 
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of a circle which cannot be parameterized by polynomial functions. Thus for applications 
where speed is more important than accuracy the Bezier curves are the favorite solution, but 
whenever accuracy is crucial or an exact solution is called for then the smooth connection 
functions are needed. 

VI. CONCLUSIONS AND OUTLOOK 

We have presented a generic solution of how to connect two points in a plane by a smooth 
curve that goes through these points with a given slope. Our approach uses extensively 
notions that are well-known in classical mechanics. The smooth connection function y(x) 
has to satisfy certain boundary conditions and to minimize the action functional that 
reflects the smoothness of any function by integrating over the inverse curvature radius (to 
some power v) times the length element along the curve. Minimizing the action via the 
Euler-Lagrange formalism leads to the equation of motion which is a complicated non-linear 
third order differential equation. However the translational and rotational symmetries of 
the problem are of great help, since they imply conserved charges, i.e. the linear and the 
angular momenta, that help to simplify the problem significantly. Making a specific choice 
for the charges it is possible to obtain the SCF y(x) for a given v in terms of hypergeometric 
functions. Applying the appropriate coordinate transformation to this solution allows to 
adjust the SCF to arbitrary boundary conditions. 

We have worked out in detail the basic formalism to find explicitly the smoothest con- 
nection between two points in the two dimensional Euclidean space. Several generalizations 
and extensions of this basic problems are possible: First, the number of points that de- 
fine the SCF can be increased. For a single curve the solution will still be determined by 
four boundary conditions, but one may choose different conditions than those in eqs. (0) 
and (0). Also it is possible to paste together several elementary solutions in order to find 
interpolations between several points which cannot be achieved by a single SCF. How to 
do this best gives rise to a new optimization problem. Finally we note that one can also 
choose to apply our formalism to a different geometry. For example one may consider a 
time-dependent SCF x(t) in Minkowski space, where the metric is defined via r 2 = x 2 — t 2 . 
Also extensions of the problem to higher dimensional spaces are conceivable. In this case 
the SCF would describe some smooth manifold that connect two extended objects (like 
strings). 
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FIGURES 




X 



FIG. 1. Different "Smooth Connection Functions" that connect the initial point P, : (xi,yi) 
with the final point Pf : (xf,yf). At Pi and Pf the curves are tangent to the associated rays 
(dashed) which are specified by their inclination angles, ai and ctf, with respect to the vector 
pointing from Pi to Pf. The light-gray curve (resulting from v = 1.2) is the shortest, but has 
the strongest curvature. Conversely, the dark-grey curve {v = 20) is the longest, but almost 
corresponds to an arc of a circle with minimal curvature. The intermediate grey curve {y = 2) 
presents a more balanced compromise between minimal length and curvature. 
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X 



FIG. 2. Standard smooth connection functions for various v. The two branches are monotonic 
and we have chosen the constants of integration such that they can be obtained from each other 
by a reflection with respect to the y-axis. Each SCF changes the sign of its curvature at x = 0. 
For large \u\ the SCF is very close to the arc of a circle. The curves for v > 1 are below the curve 
corresponding to — > oo. These curves become "flatter" and thus shorter for smaller values of 
v. They approach x = with a vanishing slope and their curvature radius diverges at x = 0. The 
closer v is to unity the sooner the curve approaches the value y(0) when x — ► 0. The curves for 
v < all reside above the curve corresponding to \v\ — > oo. For these curves both the slope q(x) 
and the curvature radius r(x) vanish at x = 0. 
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y(x) 




x 



FIG. 3. Extended standard smooth connection functions for various v. These functions are 
obtained from the SCFs of Fig. ^ by exchanging the coordinates x and y. All curves have been 
rescaled such that they are defined in the interval — 1 < x < 1. The interval —1/2 < x < 1/2 
corresponds to the SCFs of Fig. |2[ The continuations beyond x = ±1/2 use a segment of the 
other branch, which is attached to the centerpiece such that the curvature radius is continuous 
at x = ±1/2. 
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v=-2 



FIG. 4. Example for a specific set of smooth connection functions that all satisfy the same 
boundary conditions. Each curve corresponds to a different v as indicated in the legend. All the 
solutions consists of a segment of an extended standard smooth connection function that can be 
viewed as a template which may be rotated, translated and scaled in order to fit the boundary 
conditions. This example corresponds to on = 20° and aj = 70°. 



23 




v=20 v=2 v=1.3 v=-0.5 v=-2 



FIG. 5. Examples for solving specific choices of (ai,ctf) with the (extended) standard SCFs. 
Each plot corresponds to a particular choice of (ai,ctf) and shows the behavior of the solution 
for five different values of v as indicated in the legend. 
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